An ensemble penalized regression method for multi-ancestry polygenic risk prediction.

Great efforts are being made to develop advanced polygenic risk scores (PRS) to improve the prediction of complex traits and diseases. However, most existing PRS are primarily trained on European ancestry populations, limiting their transferability to non-European populations. In this article, we propose a novel method for generating multi-ancestry Polygenic Risk scOres based on enSemble of PEnalized Regression models (PROSPER). PROSPER integrates genome-wide association studies (GWAS) summary statistics from diverse populations to develop ancestry-specific PRS with improved predictive power for minority populations. The method uses a combination of L 1 (lasso) and L 2 (ridge) penalty functions, a parsimonious specification of the penalty parameters across populations, and an ensemble step to combine PRS generated across different penalty parameters. We evaluate the performance of PROSPER and other existing methods on large-scale simulated and real datasets, including those from 23andMe Inc., the Global Lipids Genetics Consortium, and All of Us. Results show that PROSPER can substantially improve multi-ancestry polygenic prediction compared to alternative methods across a wide variety of genetic architectures. In real data analyses, for example, PROSPER increased out-of-sample prediction R2 for continuous traits by an average of 70% compared to a state-of-the-art Bayesian method (PRS-CSx) in the African ancestry population. Further, PROSPER is computationally highly scalable for the analysis of large SNP contents and many diverse populations.

MUSSEL: Enhanced Bayesian polygenic risk prediction leveraging information across multiple ancestry groups.

Polygenic risk scores (PRSs) are now showing promising predictive performance on a wide variety of complex traits and diseases, but there exists a substantial performance gap across populations. We propose MUSSEL, a method for ancestry-specific polygenic prediction that borrows information in summary statistics from genome-wide association studies (GWASs) across multiple ancestry groups via Bayesian hierarchical modeling and ensemble learning. In our simulation studies and data analyses across four distinct studies, totaling 5.7 million participants with a substantial ancestral diversity, MUSSEL shows promising performance compared to alternatives. For example, MUSSEL has an average gain in prediction R2 across 11 continuous traits of 40.2% and 49.3% compared to PRS-CSx and CT-SLEB, respectively, in the African ancestry population. The best-performing method, however, varies by GWAS sample size, target ancestry, trait architecture, and linkage disequilibrium reference samples; thus, ultimately a combination of methods may be needed to generate the most robust PRSs across diverse populations.

A new method for multiancestry polygenic prediction improves performance across diverse populations.

Polygenic risk scores (PRSs) increasingly predict complex traits; however, suboptimal performance in non-European populations raise concerns about clinical applications and health inequities. We developed CT-SLEB, a powerful and scalable method to calculate PRSs, using ancestry-specific genome-wide association study summary statistics from multiancestry training samples, integrating clumping and thresholding, empirical Bayes and superlearning. We evaluated CT-SLEB and nine alternative methods with large-scale simulated genome-wide association studies (~19 million common variants) and datasets from 23andMe, Inc., the Global Lipids Genetics Consortium, All of Us and UK Biobank, involving 5.1 million individuals of diverse ancestry, with 1.18 million individuals from four non-European populations across 13 complex traits. Results demonstrated that CT-SLEB significantly improves PRS performance in non-European populations compared with simple alternatives, with comparable or superior performance to a recent, computationally intensive method. Moreover, our simulation studies offered insights into sample size requirements and SNP density effects on multiancestry risk prediction.

Effect of duck interferon-α and an anti-cap protein polyclonal antibody against duck circovirus.

Duck circovirus (DuCV) is one of the most prevalent viruses in the duck breeding industry, and causes persistent infection and severe immunosuppression. Currently, there is a serious lack of prevention and control measures and no commercial vaccine against DuCV. Therefore, effective antiviral drugs are important for treating DuCV infection. Interferon (IFN) is an important component of antiviral innate immunity, but it remains unclear whether duck IFN-α has a clinical effect against DuCV. Antibody therapy is an important way to treat viral infections. The DuCV structural protein (cap) is immunogenic, and it remains to be determined whether an anti-cap protein antibody can effectively block DuCV infection. In this study, the duck IFN-α gene and the DuCV structural protein cap gene were cloned, expressed and purified in Escherichia coli to prepare duck recombinant IFN-α and the cap protein. Then, rabbits were immunized with the recombinant cap protein to prepare a rabbit polyclonal antibody. This study investigated the antiviral effect of duck recombinant IFN-α and the anti-cap protein antibody and their combined effect on Cherry Valley ducks infected with DuCV. The results showed that the treatment significantly alleviated the clinical symptoms of immune organ atrophy and immunosuppression compared with the control. The histopathological damage of the target organs was alleviated, and replication of DuCV in the immune organs was significantly inhibited. The treatment also reduced the damage caused by DuCV to the liver and immune function, and increased the level of the DuCV antibody in the blood, thereby improving antiviral activity. Notably, the combination of duck IFN-α and the polyclonal antibody completely blocked DuCV infection after 13 days under the experimental conditions, showing a better inhibitory effect on DuCV infection than single treatments. These results showed that duck recombinant IFN-α and the anti-cap protein antibody can be used as antiviral drugs to clinically treat and control DuCV infection, particularly the vertical transmission of the virus in breeding ducks.

Hexavalent chromium disrupts the skin barrier by targeting ROS-mediated mitochondrial pathway apoptosis in keratinocytes.

Hexavalent chromium (Cr(VI)), a toxic heavy metal, is ubiquitous in daily life. Exposure to this toxic substance in occupational settings can cause dermatitis and cancer. As the body's largest organ, the skin plays a crucial role in protecting the organism against external aggressions. While previous studies have focused on the effects of Cr(VI) on skin inflammation, this study investigates the potential toxicity of Cr(VI) from the skin barrier and integrity perspective. The in vivo results of this study showed that mice exposed to Cr(VI) experienced skin deterioration and hemorrhaging, as well as a reduction in the thickness of the collagen fiber layer. TUNEL and Occludin staining results revealed that Cr(VI)'s toxicity primarily targeted keratinocytes. Experiments in vitro demonstrated that Cr(VI) treatment decreased the activity of HaCaT cells, altered cell morphology, and increased LDH secretion. Further research revealed that Cr(VI) could modify membrane permeability, impair membrane integrity, and reduce the protein expression of ZO-1 and Occludin. In addition, it was discovered that Cr(VI) promoted cell apoptosis and inhibited AKT activation. However, the addition of a caspase inhibitor and an AKT activator prevented Cr(VI)-induced injury to the cell membrane barrier, indicating that apoptosis plays a crucial role in this process. The addition of three apoptotic pathway inhibitors, confirmed that Cr(VI) damaged the cell barrier through ROS-mediated mitochondrial pathway apoptosis. Moreover, the use of a ROS inhibitor significantly reduced Cr(VI)-induced apoptosis and cell barrier injury. In conclusion, this study provides an experimental foundation for the treatment of skin injury caused by Cr(VI).

Comparison study of protective effects of porcine bile acids and sheep bile acids against heat stress in chickens.

BACKGROUND: Heat stress (HS) is known to exert negative effects on the poultry and breeding industry, resulting in severe economic losses. Bile acids (BAs), an important component of bile, play a crucial role in improving the production performance of livestock and poultry, alleviating stress injury, and ensuring the health of livestock and poultry. At present, porcine BAs are widely used because of their therapeutic effects on HS; however, it remains unclear whether the same effects are exerted by sheep BAs, which are different from porcine BAs and have different compositions. In this study, we compared the anti-HS effects of porcine BAs and sheep BAs in the diet by establishing an HS model of chicks and investigating the chicken performance, HS-related genes' expression, oxidative stress markers, jejunal histoarchitecture, inflammatory cytokines' expression, jejunal secreted immunoglobulin A concentration, and cecal bacterial flora. RESULTS: The results showed that the addition of sheep BAs to the diet increased the average daily weight gain and the feed conversion ratio of chicks. Under HS, sheep BAs were more effective than porcine BAs in improving the activities of lactate dehydrogenase and glutamic pyruvic transaminase in serum and the content/activity of malondialdehyde, superoxide dismutase, and reduced glutathione in serum and tissue, in reducing the messenger RNA (mRNA) expression of heat shock proteins (HSP60, HSP70, and HSP90) in the liver and jejunum, and in improving the histological structure and the expression of tight junction proteins (occludin and zonula occludens-1) and enriching intestinal bacterial flora. However, porcine BAs were significantly inferior to sheep BAs in reducing the mRNA expression of inflammatory factors (interleukin-6, interleukin-1ß, and tumor necrosis factor-α). CONCLUSION: The effect of sheep BAs was more significant than porcine BAs was in alleviating HS injury in chicks, suggesting that sheep BAs have great potential as new feed nutrition and health additive to improve poultry production performance and prevent HS. © 2023 Society of Chemical Industry.

MUSSEL: Enhanced Bayesian Polygenic Risk Prediction Leveraging Information across Multiple Ancestry Groups.

Polygenic risk scores (PRS) are now showing promising predictive performance on a wide variety of complex traits and diseases, but there exists a substantial performance gap across different populations. We propose MUSSEL, a method for ancestry-specific polygenic prediction that borrows information in the summary statistics from genome-wide association studies (GWAS) across multiple ancestry groups. MUSSEL conducts Bayesian hierarchical modeling under a MUltivariate Spike-and-Slab model for effect-size distribution and incorporates an Ensemble Learning step using super learner to combine information across different tuning parameter settings and ancestry groups. In our simulation studies and data analyses of 16 traits across four distinct studies, totaling 5.7 million participants with a substantial ancestral diversity, MUSSEL shows promising performance compared to alternatives. The method, for example, has an average gain in prediction R2 across 11 continuous traits of 40.2% and 49.3% compared to PRS-CSx and CT-SLEB, respectively, in the African Ancestry population. The best-performing method, however, varies by GWAS sample size, target ancestry, underlying trait architecture, and the choice of reference samples for LD estimation, and thus ultimately, a combination of methods may be needed to generate the most robust PRS across diverse populations.

An Ensemble Penalized Regression Method for Multi-ancestry Polygenic Risk Prediction.

Great efforts are being made to develop advanced polygenic risk scores (PRS) to improve the prediction of complex traits and diseases. However, most existing PRS are primarily trained on European ancestry populations, limiting their transferability to non-European populations. In this article, we propose a novel method for generating multi-ancestry Polygenic Risk scOres based on enSemble of PEnalized Regression models (PROSPER). PROSPER integrates genome-wide association studies (GWAS) summary statistics from diverse populations to develop ancestry-specific PRS with improved predictive power for minority populations. The method uses a combination of ℒ1 (lasso) and ℒ2 (ridge) penalty functions, a parsimonious specification of the penalty parameters across populations, and an ensemble step to combine PRS generated across different penalty parameters. We evaluate the performance of PROSPER and other existing methods on large-scale simulated and real datasets, including those from 23andMe Inc., the Global Lipids Genetics Consortium, and All of Us. Results show that PROSPER can substantially improve multi-ancestry polygenic prediction compared to alternative methods across a wide variety of genetic architectures. In real data analyses, for example, PROSPER increased out-of-sample prediction R2 for continuous traits by an average of 70% compared to a state-of-the-art Bayesian method (PRS-CSx) in the African ancestry population. Further, PROSPER is computationally highly scalable for the analysis of large SNP contents and many diverse populations.

Pathogenicity of duck circovirus and fowl adenovirus serotype 4 co-infection in Cherry Valley ducks.

Duck circovirus (DuCV) is one of the most prevalent infectious viruses in the duck industry in China. Although the clinical symptoms vary, it often causes immunosuppression in the host and leads to secondary infection with other pathogens. Fowl adenovirus serotype 4 (FAdV-4) mainly infects chickens and causes hydropericardium hepatitis syndrome. However, the incidence of infection in ducks has increased in recent years, and the phenomenon of mixed infection with DuCV is very common, resulting in more severe clinical morbidity. However, there is no systematic study evaluating the presence of mixed infection. To explore the synergistic pathogenicity of DuCV and FAdV-4 co-infection in Cherry Valley ducks, a comparative experiment was established between DuCV and FAdV-4 co-infection and single infection animal models. It was found that DuCV and FAdV-4 co-infected ducks showed more pronounced clinical signs of pericardial effusion, hepatitis and immunosuppression; more severe tissue damage in target organs; and more significant levels of viral load, biochemical indicators and immune indicators in various organs compared with Cherry Valley ducks infected with just one virus. The results showed that co-infection with DuCV and FAdV-4 may promote greater viral replication, causing more severe tissue damage and immunosuppression than infection with just one virus. Therefore, the monitoring and prevention of the two viruses should be strengthened clinically, with a particular focus on the potential harm of DuCV as it carries the highest infection rate.

A saturated map of common genetic variants associated with human height.

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

Evaluation of the antiviral effect of four plant polysaccharides against duck circovirus.

Recently, outbreaks of duck circovirus (DuCV) are frequently occurring worldwide due to secondary infections caused by post infection-induced immunosuppression. Due to a lack of preventive drugs and vaccines, the waterfowl industry losses are ever increasing. In this study, we extracted Astragalus polysaccharides (APS), pine pollen polysaccharides (PPPS), Aloe vera polysaccharides (AVE), and Ficus carica polysaccharides (FCPS) from Astragalus, pine pollen, aloe, and F. carica leaves, respectively. We randomly divided 150 one-day-old Cherry Valley ducks into five groups, which were inoculated with the DuCV solution and orally administered APS, PPPS, AVE, FCPS, and phosphate buffer saline (PBS), respectively. We collected the duck immune organs and serum samples at 8, 16, 24, 32, 40, and 48 days post-infection (dpi). Using clinical symptom analysis, molecular biology experiments, and serological experiments, we proved that plant polysaccharides could (a) improve the duck immunity, (b) reduce the viral load, and (c) mitigate DuCV-induced damage to immune organs, with both APS and PPPS having significant effects. Moreover, we detected viral load and cytokines within the first 8 dpi. Since the body's innate immunity could inhibit viral replication within five days of virus infection, 1-5 dpi was the best treatment time. Among the four polysaccharides showing in vitro anti-apoptotic activity, APS and PPPS significantly inhibited the DuCV infection-induced apoptosis of peripheral blood lymphocytes. Overall, since our findings show APS and PPPS having significant anti-DuCV effects both in vivo and in vitro, they can be promising candidates for preventing DuCV infection in ducks.

Polygenic prediction of educational attainment within and between families from genome-wide association analyses in 3 million individuals.

We conduct a genome-wide association study (GWAS) of educational attainment (EA) in a sample of ~3 million individuals and identify 3,952 approximately uncorrelated genome-wide-significant single-nucleotide polymorphisms (SNPs). A genome-wide polygenic predictor, or polygenic index (PGI), explains 12-16% of EA variance and contributes to risk prediction for ten diseases. Direct effects (i.e., controlling for parental PGIs) explain roughly half the PGI's magnitude of association with EA and other phenotypes. The correlation between mate-pair PGIs is far too large to be consistent with phenotypic assortment alone, implying additional assortment on PGI-associated factors. In an additional GWAS of dominance deviations from the additive model, we identify no genome-wide-significant SNPs, and a separate X-chromosome additive GWAS identifies 57.

Characterizing mood disorders in the AFFECT study: a large, longitudinal, and phenotypically rich genetic cohort in the US.

There has recently been marked progress in identifying genetic risk factors for major depression (MD) and bipolar disorder (BD); however, few systematic efforts have been made to elucidate heterogeneity that exists within and across these diagnostic taxa. The Affective disorders, Environment, and Cognitive Trait (AFFECT) study presents an opportunity to identify and associate the structure of cognition and symptom-level domains across the mood disorder spectrum in a prospective study from a diverse US population.Participants were recruited from the 23andMe, Inc research participant database and through social media; self-reported diagnosis of MD or BD by a medical professional and medication status data were used to enrich for mood-disorder cases. Remote assessments were used to acquire an extensive range of phenotypes, including mood state, transdiagnostic symptom severity, task-based measures of cognition, environmental exposures, personality traits. In this paper we describe the study design, and the demographic and clinical characteristics of the cohort. In addition we report genetic ancestry, SNP heritability, and genetic correlations with other large cohorts of mood disorders.A total of 48,467 participants were enrolled: 14,768 with MD, 9864 with BD, and 23,835 controls. Upon enrollment, 47% of participants with MD and 27% with BD indicated being in an active mood episode. Cases reported early ages of onset (mean = 13.2 and 14.3 years for MD and BD, respectively), and high levels of recurrence (78.6% and 84.9% with >5 episodes), psychotherapy, and psychotropic medication use. SNP heritability on the liability scale for the ascertained MD participants (0.19-0.21) was consistent with the high level of disease severity in this cohort, while BD heritability estimates (0.16-0.22) were comparable to reports in other large scale genomic studies of mood disorders. Genetic correlations between the AFFECT cohort and other large-scale cohorts were high for MD but not for BD. By incorporating transdiagnostic symptom assessments, repeated measures, and genomic data, the AFFECT study represents a unique resource for dissecting the structure of mood disorders across multiple levels of analysis. In addition, the fully remote nature of the study provides valuable insights for future virtual and decentralized clinical trials within mood disorders.

Multi-Trait Genetic Analysis Identifies Autoimmune Loci Associated with Cutaneous Melanoma.

Genome-wide association studies (GWAS) have identified a number of risk loci for cutaneous melanoma. Cutaneous melanoma shares overlapping genetic risk (genetic correlation) with a number of other traits, including its risk factors such as sunburn propensity. This genetic correlation can be exploited to identify additional cutaneous melanoma risk loci by multitrait analysis of GWAS (MTAG). We used bivariate linkage disequilibrium-score regression score regression to identify traits that are genetically correlated with clinically confirmed cutaneous melanoma and then used publicly available GWAS for these traits in a multitrait analysis of GWAS. Multitrait analysis of GWAS allows GWAS to be combined while accounting for sample overlap and incomplete genetic correlation. We identified a total of 74 genome-wide independent loci, 19 of them were not previously reported in the input cutaneous melanoma GWAS meta-analysis. Of these loci, 55 were replicated (P < 0.05/74, Bonferroni-corrected P-value in two independent cutaneous melanoma replication cohorts from Melanoma Institute Australia and 23andMe, Inc. Among the, to our knowledge, previously unreported cutaneous melanoma loci are ones that have also been associated with autoimmune traits including rs715199 near LPP and rs10858023 near AP4B1. Our analysis indicates genetic correlation between traits can be leveraged to identify new risk genes for cutaneous melanoma.

The Genetic Architecture of Depression in Individuals of East Asian Ancestry: A Genome-Wide Association Study.

Importance: Most previous genome-wide association studies (GWAS) of depression have used data from individuals of European descent. This limits the understanding of the underlying biology of depression and raises questions about the transferability of findings between populations. Objective: To investigate the genetics of depression among individuals of East Asian and European descent living in different geographic locations, and with different outcome definitions for depression. Design, Setting, and Participants: Genome-wide association analyses followed by meta-analysis, which included data from 9 cohort and case-control data sets comprising individuals with depression and control individuals of East Asian descent. This study was conducted between January 2019 and May 2021. Exposures: Associations of genetic variants with depression risk were assessed using generalized linear mixed models and logistic regression. The results were combined across studies using fixed-effects meta-analyses. These were subsequently also meta-analyzed with the largest published GWAS for depression among individuals of European descent. Additional meta-analyses were carried out separately by outcome definition (clinical depression vs symptom-based depression) and region (East Asian countries vs Western countries) for East Asian ancestry cohorts. Main Outcomes and Measures: Depression status was defined based on health records and self-report questionnaires. Results: There were a total of 194 548 study participants (approximate mean age, 51.3 years; 62.8% women). Participants included 15 771 individuals with depression and 178 777 control individuals of East Asian descent. Five novel associations were identified, including 1 in the meta-analysis for broad depression among those of East Asian descent: rs4656484 (ß = -0.018, SE = 0.003, P = 4.43x10-8) at 1q24.1. Another locus at 7p21.2 was associated in a meta-analysis restricted to geographically East Asian studies (ß = 0.028, SE = 0.005, P = 6.48x10-9 for rs10240457). The lead variants of these 2 novel loci were not associated with depression risk in European ancestry cohorts (ß = -0.003, SE = 0.005, P = .53 for rs4656484 and ß = -0.005, SE = 0.004, P = .28 for rs10240457). Only 11% of depression loci previously identified in individuals of European descent reached nominal significance levels in the individuals of East Asian descent. The transancestry genetic correlation between cohorts of East Asian and European descent for clinical depression was r = 0.413 (SE = 0.159). Clinical depression risk was negatively genetically correlated with body mass index in individuals of East Asian descent (r = -0.212, SE = 0.084), contrary to findings for individuals of European descent. Conclusions and Relevance: These results support caution against generalizing findings about depression risk factors across populations and highlight the need to increase the ancestral and geographic diversity of samples with consistent phenotyping.

FUT6 deficiency compromises basophil function by selectively abrogating their sialyl-Lewis x expression.

Sialyl-Lewis x (sLex, CD15s) is a tetra-saccharide on the surface of leukocytes required for E-selectin-mediated rolling, a prerequisite for leukocytes to migrate out of the blood vessels. Here we show using flow cytometry that sLex expression on basophils and mast cell progenitors depends on fucosyltransferase 6 (FUT6). Using genetic association data analysis and qPCR, the cell type-specific defect was associated with single nucleotide polymorphisms (SNPs) in the FUT6 gene region (tagged by rs17855739 and rs778798), affecting coding sequence and/or expression level of the mRNA. Heterozygous individuals with one functional FUT6 gene harbor a mixed population of sLex+ and sLex- basophils, a phenomenon caused by random monoallelic expression (RME). Microfluidic assay demonstrated FUT6-deficient basophils rolling on E-selectin is severely impaired. FUT6 null alleles carriers exhibit elevated blood basophil counts and a reduced itch sensitivity against insect bites. FUT6-deficiency thus dampens the basophil-mediated allergic response in the periphery, evident also in lower IgE titers and reduced eosinophil counts.

Genetic predisposition to mosaic Y chromosome loss in blood.

Mosaic loss of chromosome Y (LOY) in circulating white blood cells is the most common form of clonal mosaicism1-5, yet our knowledge of the causes and consequences of this is limited. Here, using a computational approach, we estimate that 20% of the male population represented in the UK Biobank study (n = 205,011) has detectable LOY. We identify 156 autosomal genetic determinants of LOY, which we replicate in 757,114 men of European and Japanese ancestry. These loci highlight genes that are involved in cell-cycle regulation and cancer susceptibility, as well as somatic drivers of tumour growth and targets of cancer therapy. We demonstrate that genetic susceptibility to LOY is associated with non-haematological effects on health in both men and women, which supports the hypothesis that clonal haematopoiesis is a biomarker of genomic instability in other tissues. Single-cell RNA sequencing identifies dysregulated expression of autosomal genes in leukocytes with LOY and provides insights into why clonal expansion of these cells may occur. Collectively, these data highlight the value of studying clonal mosaicism to uncover fundamental mechanisms that underlie cancer and other ageing-related diseases.

Genetic Architectures of Childhood- and Adult-Onset Asthma Are Partly Distinct.

The extent to which genetic risk factors are shared between childhood-onset (COA) and adult-onset (AOA) asthma has not been estimated. On the basis of data from the UK Biobank study (n = 447,628), we found that the variance in disease liability explained by common variants is higher for COA (onset at ages between 0 and 19 years; h2g = 25.6%) than for AOA (onset at ages between 20 and 60 years; h2g = 10.6%). The genetic correlation (rg) between COA and AOA was 0.67. Variation in age of onset among COA-affected individuals had a low heritability (h2g = 5%), which we confirmed in independent studies and also among AOA-affected individuals. To identify subtype-specific genetic associations, we performed a genome-wide association study (GWAS) in the UK Biobank for COA (13,962 affected individuals) and a separate GWAS for AOA (26,582 affected individuals) by using a common set of 300,671 controls for both studies. We identified 123 independent associations for COA and 56 for AOA (37 overlapped); of these, 98 and 34, respectively, were reproducible in an independent study (n = 262,767). Collectively, 28 associations were not previously reported. For 96 COA-associated variants, including five variants that represent COA-specific risk factors, the risk allele was more common in COA- than in AOA-affected individuals. Conversely, we identified three variants that are stronger risk factors for AOA. Variants associated with obesity and smoking had a stronger contribution to the risk of AOA than to the risk of COA. Lastly, we identified 109 likely target genes of the associated variants, primarily on the basis of correlated expression quantitative trait loci (up to n = 31,684). GWAS informed by age of onset can identify subtype-specific risk variants, which can help us understand differences in pathophysiology between COA and AOA and so can be informative for drug development.

Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero.

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sex-specific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N = 15 661, with replication N = 75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (ß = 0.06; P = 0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.

Powerful and robust cross-phenotype association test for case-parent trios.

There has been increasing interest in identifying genes within the human genome that influence multiple diverse phenotypes. In the presence of pleiotropy, joint testing of these phenotypes is not only biologically meaningful but also statistically more powerful than univariate analysis of each separate phenotype accounting for multiple testing. Although many cross-phenotype association tests exist, the majority of such methods assume samples composed of unrelated subjects and therefore are not applicable to family-based designs, including the valuable case-parent trio design. In this paper, we describe a robust gene-based association test of multiple phenotypes collected in a case-parent trio study. Our method is based on the kernel distance covariance (KDC) method, where we first construct a similarity matrix for multiple phenotypes and a similarity matrix for genetic variants in a gene; we then test the dependency between the two similarity matrices. The method is applicable to either common variants or rare variants in a gene, and resulting tests from the method are by design robust to confounding due to population stratification. We evaluated our method through simulation studies and observed that the method is substantially more powerful than standard univariate testing of each separate phenotype. We also applied our method to phenotypic and genotypic data collected in case-parent trios as part of the Genetics of Kidneys in Diabetes (GoKinD) study and identified a genome-wide significant gene demonstrating cross-phenotype effects that was not identified using standard univariate approaches.